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Abstract 

In charged current weak interaction processes, neutrinos are produced in an entangled state with 

^"•^ the charged lepton. This correlated state is disentangled by the measurement of the charged lepton 

O ill ^ detector at the production site. We study the dynamical aspects of disentanglement, propa- 

^^ gation and detection, in particular the conditions under which the disentangled state is a coherent 

^ superposition of mass eigenstates. The appearance and disappearance far-detection processes are 

described from the time evolution of this disentangled "collapsed" state. The familiar quantum 

mechanical interpretation and factorization of the detection rate emerges when the quantum state 

is disentangled on time scales much shorter than the inverse oscillation frequency, in which case 

the final detection rate factorizes in terms of the usual quantum mechanical transition probability 

Q^ provided the final density of states is insensitive to the neutrino energy difference. We suggest 

Q^ possible corrections for short-baseline experiments. If the charged lepton is unobserved, neutrino 

V oscillations and coherence are described in terms of a reduced density matrix obtained by tracing 

I— I out an un-observed charged lepton. The diagonal elements in the mass basis describe the pro- 

P^ duction of mass eigenstates whereas the off diagonal ones provide a measure of coherence. It is 

J> shown that coherences are of the same order of the diagonal terms on time scales up to the inverse 

S:^ oscillation frequency, beyond which the coherences oscillate as a result of the interference between 

^ mass eigenstates. 
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I. INTRODUCTION 

Neutrino masses, mixing and oscillations are the clearest evidence yet of physics beyond 
the standard model [TH3]. They provide an explanation of the solar neutrino problem |3]- 
El and have other important phenomenological [IHSl rrHTP]- astrophysical [TT], [T2] and 
cosmological [13] consequences. Furthermore, they offer an incredible example of macro- 
scopic quantum coherence, where in long baseline oscillation experiments this coherence is 
maintained over hundreds of kilometers. 

Given the importance of neutrino oscillations in pointing us to the next layer beyond 
that of the standard model, it is not unreasonable to make sure we understand the domain 
of validity of the various calculations of oscillation probabilities. In particular, the standard 
approach of treating neutrino oscillations by analogy with Rabi-oscillations in a two state 
system (see for example, P-El El E] and references therein) while simple and intuitive, is not 
necessarily the last word on the subject. Deeper investigations of this basic paradigm have 
already raised a number of important and fundamental questions [HHlT] that are still being 
debated ^Ml\. 

If we truly want to understand the results from oscillation experiments, we need to 
factor in both the production and detection mechanisms. In term of its detection, the 
neutrino is never directly detected. We only know that it is there due to the charged leptons 
associated with its weak interaction at the detector. As for the production mechanism, the 
neutrino state is produced by the decay of a parent particle via charged current interactions. 
Interestingly, a discussion of coherence aspects of neutrinos [22H2Z| has recognized that the 
neutrino state produced by the decay of a parent particle via charged current interactions 
is in fact entangled with that of the charged lepton. 

How does this entangled state evolve? As usual, there are two choices. If we measure 
the charged lepton, this disentangles the quantum state at a time determined either by the 
detection process of the charged lepton at the source or by stopping at a "beam dump" near 
the source. This disentangled neutrino state then evolves further in time and is detected at 
the "far detector" via another charged lepton. Therefore the production, disentanglement 
and detection involve two times: the time at which the charged lepton is observed at the 
source (if observed at all) or stopped at a beam dump near the source, which is when the 
neutrino/charged lepton state is disentangled, and the time when the neutrino is observed 
via a charged lepton in the far detector. We note that for long-baseline experiments these 
two time scales are widely separated. 

If the charged lepton is not observed, the state is described by a reduced density matrix, 
which is obtained by tracing out the charged lepton. In this case, we can track the behavior of 
the off-diagonal elements, i.e., the coherences which contain information on the interference 
between the different mass eigenstates. 

What we do in this article is to use a full field theoretic approach to follow this time 
evolution explicitly. We start from a simplified model that encodes the salient interactions 
necessary to analyze the development of oscillations. Using this model, we calculate the 
form of the state produced when the parent particle (be it a. W boson or a pion) decays 
into the flavor neutrino and its associated charged lepton. We then use the time evolution 
operator, calculated within perturbation theory from the interaction Lagrangian, to evolve 
this state until it reaches the detector. At this point, we obtain the overlap with the final 
states observed at the detector. 

Our main observation is that oscillation experiments are two-time experiments, the two 



times being the time at which the original charged lepton is measured and that when the final 
one is measured. While this does not seem earth-shattering, what is extremely interesting 
is the interplay between these times (and/or time differences) and other time scales in the 
problem such as the oscillation time scale and the time scale upon which it can be said that 
the two mass eigenstates have been resolved. In fact, one of our results is that the treatment 
of the flavor eigenstate as the usual coherent mixture of the two mass eigenstates requires 
that the original entangled state already be disentangled long before an oscillation time. 
While this is true for long baseline experiments, it is not necessarily so for short baseline 
ones, and this could give rise to a reinterpretation of their results to take the dependence 
on the disentangling time into account. 

In all of this, we would like to emphasize that we are just following the tenets of quantum 
mechanics in our calculation: prepare a state, construct and then use the relevant time 
evolution operator and lastly, make a measurement on that state by overlapping it with 
various states of interest at the later time. Following the time evolution in detail as we do 
gives us the access to flner details of the behavior of the state than what would be obtained 
when one uses Fermi's Golden Rule, which requires taking an inflnite time limit, to directly 
calculate rates and the like. In our formalism, the flnite-time effects become explicit and this 
allows us to fully understand the domain of validity of the various results in the literature. 



II. A MODEL OF "NEUTRINO" OSCILLATIONS 

Let us start by introducing our model. In order to exhibit the main results in a clear and 
simple manner, we introduce a bosonic model that describes mixing, oscillations and charged 
current weak interactions reliably without the complications associated with fermionic and 
gauge fields. We can do so because the technical complications associated with spinors and 
gauge fields are irrelevant to the physics of mixing and oscillations, as is obviously manifest 
in meson mixing. Our model is defined by the following Lagrangian density 

£ = £o[W^,U + 'CoH + £int[W^,L,'^a] ; a = e,/i (2.1) 

with 

Cq[u] = ^ [9^^^9^^ - ^"^M^j , (2.2) 

where \1/ is a flavor doublet representing the neutrinos 

^ = f !^0 , (2.3) 

and M is the mass matrix 



^M 
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(2.4) 



The interaction Lagrangian is similar to the charged current interaction of the standard 
model, namely 

Cint{x,t) = gW{x,t) le{x,t)ue{x,t) + lf,{x,t)Uf,{x,t) , (2.5) 

where g is the coupling constant. W{x) represents the vector boson, or alternatively the 
pion fleld, and /„, a = e,fi the two charged leptons. The mass matrix is diagonalized by a 



unitary transformation 

U'\e) M ,m = ( 7 I ) ; U(0) = ^ :~", — j , (2,6) 

In terms of the doublet of mass eigenstates, the flavor doublet can be expressed as 
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(2.7) 



This bosonic model clearly describes charged current weak interactions reliably as it includes 
all the relevant aspects of mixing and oscillations. 

We consider "neutrino" oscillation experiments following the interaction processes illus- 
trated in Fig. [T| 
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FIG. 1: Typical experiment in which the charged leptons are measured at the source and far 
detector and the neutrino is an intermediate state. 



Throughout this work, we will use plane wave states for simplicity of exposition, though 
we discuss some of the modifications that the use of wave packets would require later in this 
section. 

In order to study aspects of coherence we consider a simplified interaction Lagrangian 
density 

Ci = gWeVe = 5'PFe(cos6'z/i + sin6'z/2), (2.9) 

focusing only on one lepton, which we refer to as the "electron." The full coupling as in Eq. 



(2.5) can be treated similarly without modifying the main conclusions. Although W may 



be interpreted as a charged vector boson, the analysis is obviously the same if it describes a 
pion field. 

We can study aspects of coherence by focusing on the Fock state obtained upon evolution 
of the decaying initial state. We consider a plane-wave Fock initial state |iy(A;)) at ti = 0. 
The time evolution operator is 



^-iH{t-U) ^ ^-iHot 



U{t,U 



JHnti 



(2.10) 



where 



Uit,ti 



T e 



i J^ (1? X dt Ci^t{x ,t) 



(2.11) 



is the time evolution operator in the interaction picture and Hq is the non-interacting Hamil- 
tonian. 

By time t, the initial state has evolved into |^(t)) = |iy(A;)) e~* fc + |\l'e(t)). To lowest 
order in the interaction, we find the second term to be 



\^ ^{t)) = ig e-'"°' / dt' £x W{x,t')e{x,t'){cos9iyi{x,t') + smeiy2ix,t')) \W{k)), 

(2.12) 
where all the fields are in the interaction picture. Though the field operator W can either 
annihilate the initial state or create another W particle, the state with two W particles 
features faster oscillations that will average out. In what follows, we consider only the Fock 
state resulting from the annihilation of the initial particle, leading to the state 
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in which the electron and the neutrinos are entangled} . The neutrino state that is entangled 



with the muon is obtained from (2.13) by replacing cos 6' — )■ — sin^ ; sin^ — )■ cos^. 



In what follows, we consider ultrarelativistic neutrinos and write 



VLi=VL- ^ ; ^2 = fi + A, 



(2.14) 



where 
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taking A ^ f] as is the case for ultrarelativistic nearly degenerate neutrinos. 



(2.15) 



^ The result for the wavefunction in Ref. |22j may be understood using a (non-perturbative) Wigner- 
Weisskopf approximation for the decaying parent particle, replacing Ey/ — >■ Ey/ — iVw- Taking t ^ 1/Tw 
in the integral replaces the brackets in (2.131 by \/{Ew — E^ — flj — iTw) whose absolute value is 
proportional to S{Ew — £"^ — flj)/Tw 



A. Unobserved daughter particles: time evolution of the density matrix 



If the electrons (or daughter particle in Ref. [22]) are not observed, they can be traced 
out of the density matrix |\l/e(t))(\l/e(^)|- This gives the reduced density matrix 



Pr.{t) = Tre|^,(t))(^e(t)| 
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(2.16) 



This expression contains remarkable information. The function sin^(xt)/x^ is the usual 
"diffraction" function of Fermi's Golden rule. In the formal long time limit sin^(xt)/x^ — )■ 
TTt6{x), the first two terms of the density matrix, which are the diagonal entries in the mass 
basis, describe the production process of the mass eigenstates. As will be seen below, in the 
long time limit, the time derivative of these two terms yields the production rate for the 
mass eigenstates. In the formal t — )■ oo limit, these are the diagonal terms obtained in Ref. 
[22] . where in that reference, the product of delta functions is again understood as the total 
time elapsed times an energy conserving delta function. 



The off-diagonal terms in the last line of (2.16) describe the "coherences" and display 
the oscillatory phases from the interference of the mass eigenstates. The time dependent 
factors of the off-diagonal density matrix elements determine the coherence between the mass 
eigenstates and will be a ubiquitous contribution in the real time description of oscillations 
that follows below. The functions 
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(2.17) 



are strongly peaked at x ± A = with height t and width ~ 27r/t. In the infinite time limit, 
/±(x,t. A) — )■ 211 5 {x ± A) and therefore their product would vanish in this limit, leading 
to the vanishing of the coherence. This is the result obtained in Ref. [22]. However, at 
finite time t, they feature a non-vanishing overlap when 2A < 27r/t. We recognize this as 
the condition for oscillations. We note that t ~ vr/A yields a macroscopically large time 
scale. The functions f±{x,t,A) and their products are depicted in Figs. (2][3) for the values 
A = 0.1, t = 40, 100, respectively. 



-10 
-0.6 





2A 


— 


/'\ r\ 


: f.^(x,40,0.1)/ 

1 

- 1 

1 

t 


\ \f_(x,40,0.1) 


1 

— 1 

1 


w 1 \/\ vy,^'' 






-0.4 -0.2 



0.0 
X 



0.2 0.4 




FIG. 2: The functions f±{x, t, A) vs. x for t = 40, 100, A = 0.1 
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FIG. 3: The products f-{x,t,A)f+{x,t,A) vs. x ioi t = 40, 100,A = 0.1 



It is straightforward to find 
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sin(At) 




(2.18) 



In the long time limit, the terms in the first bracket yield a sum of delta functions at x = ±A. 
Upon integrating the product of /-/+ with functions of compact support, the contribution 



from the second line in (2.18) is negligible in the long time limit. Therefore, the long time 



limit of /+/- can be replaced by 



f_{x,t,A)U{x,t,A) = 7r 



sin(At) 

A 



6{x - A) + 6{x + A) 



(2.19) 



During a time interval t ^ 27r/A, the product f-{x,t, A)/+(a;, t. A) ~ 7Tt[6{x—A)+6{x+A)] 
grows linearly in time until it begins to oscillate with frequency 27r/A for t > 27r/A. 



Therefore, we conclude that upon integration with a smooth density of states, the off- 
diagonal terms in the density matrix grow linearly in time for t <^ tosc = 27r/A, but feature 
a bound oscillatory behavior of frequency A for t > 27r/A. 



The diagonal terms, i.e., the first two terms in the reduced density matrix (2.16), are 



proportional to 4sin^(x|)/x^ — )■ 27r t5{x). This shows that the coherences or off diagonal 
terms are linear in time and of the same order as the diagonal elements for t < tosc = 27r/A, 
but are of 0{1/At) and oscillate compared to the diagonal terms for t ^ tosc- This behavior 
is displayed in Fig. (Il]), where as an example we consider a smooth density of states and 
the integral 
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The case A = describes either of the diagonal terms and is linearly secular in time. This 
figure clearly shows the slow oscillations for t > 1/A. 
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FIG. 4: The integral I{t,A) = J^^e"''^ f+{x,t,A)f^{x,t,A)dx vs. t for A = 0,0.1,0.3 



Therefore, the approximation (2.19) is reliable in the long time limit and upon in- 



tegration with functions of compact support. We see that for large t, but tA ^ 1, 

/+(x,t,A)/_(x,t,A) -^ Tit \6{x - A) + 6{x + A)] and for A -^0, the product yields 27rt(5(x). 

The reduced density matrix (2.16) allows us to obtain the time evolution of the neutrino 



populations and coherences, namely 

rii {p, t) = Tr pr {t)al (p) ai (p) ; dj {p, t) = Tr p^ (t) aj (p) a^ {p) , i^j 



(2.21) 



where the annihilation and creation operators are in the Schroedinger picture. In the long 
time limit and using the results above, we find 



^i(P) t) = t Fi(p) cos^ 9 ; n2{p, t) = t F2(p) sin^ 9, 



where 
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are the partial widths for the decay of the W into the charged lepton and the neutrino mass 
eigenstates. For the off-diagonal matrix elements or coherences we find, 
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dni{p, t)/dt yields the production rate of the neutrino mass eigenstates from the decay of 
the W and the coherences Cij are non- vanishing at any finite time. In Ref. [22], these coher- 
ences vanish as a consequence of the product of delta functions on the different mass shells of 
the mass eigenstates. The coherences vanish in the formal infinite time limit because of the 
oscillatory behavior averages out on time scales t ^ 1/A and they are of (9(sin(tA)/tA) 
with respect to the diagonal terms. However, we emphasize that experimentally, the time 
scales involved (or length scales) are of order 1/A as these are the scales on which oscillatory 
phenomena are revealed. Taking fii ~ ^^2 ~ fi in the denominators in ( 2.23[2.24 ), it follows 

that 

^ ,^ , sin20sinftAl ,a+ 1 r i / n 

Ci2(p,t) ^ — ^ e^^* 2 [^^ + ^']- ^^-^^^ 

Therefore, the coherences are of the same order of the population terms on time scales 
t < 1/A, but average out for t ^ 1/A, showing that coherence persists over the oscillation 
time scale. 



B. Disentangling the neutrino: a two-time measurement 



As we have discussed above, a long baseline experiment is actually a two-time measure- 
ment process, as the charged lepton produced at the interaction vertex at the source is 
detected at the source or stopped and absorbed in a nearby "beam dump" . 

This "measurement" of the charged lepton disentangles the neutrinos from the charged 



-iE^-^tc 



lepton in the quantum state (2.13) [22|. The detection of the charged lepton at the source (or 
its stopping at a nearby "beam dump") projects the quantum state (2.13) at the observation 

time ts onto the single particle charged lepton state e 

states into the "collapsed" state 
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We note that up to a phase the coefficient functions that multiply [1^1,2) are 
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respectively. In the limit ts — )■ 00, these become 27r6{Es—^i) and 27t6{Es—^2), respectively. 
Therefore, in this limit, for a fixed Es, one of the quantum states will be projected out. 
However, as we insist on keeping a finite time interval, we will keep ts finite. 
The state \Ve{Q,ts)) then evolves forward in time with the full Hamiltonian 



The "free" evolution is obtained by setting to lowest order U{t,ts) = 1, leading to 



(2.28) 
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The phase factors e"*^-'*'^'^ multiplying each mass eigenstate are the consequence of the 
phase build-up during the time evolution from the production vertex until the detection of 
the charged lepton and the collapse of the wave function. These can be absorbed into the 
definition of the states |i^i,2) at the fixed time ts- 
The expression (2.29) features the factors 
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which as ^5 — 7- 00 becomes 2tt 6{Es — ^j). These factors, which are a direct consequence of 
the neutrino state being produced by the decay of the "parent" particle (here the W) into 



an entangled quantum state, distinguish Eq. (2.29) from the familiar quantum mechanical 



description. These factors emerge from the (approximate) energy conservation at the decay 
vertex. Again, in the t^ — )■ 00 limit, if the energy of the parent particle and the charged 
lepton are both certain, only one of the mass eigenstates will be produced but not both. 
However, writing f2i^2 as in Eq. (2.14), it follows that for ts ^ 4Q/6m'^, the width of 



the "diffraction" functions is much larger than the frequency difference A and there is a 
substantial overlap between these "approximate" delta functions (see Fig. (|2])). Only for 
ts > tosc = 27r/A are the two peaks at Es — ^ = ±A actually resolved, whereas for ts ^ tosc, 
the two peaks are unresolved, "blurred" into one broad peak at fl. Thus for ts <^ 27r/A, we 
can use the approximation 
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To illustrate the validity of the above approximation, let us consider the case in which 
the typical size of the source or beam dump is a few meters across. In a typical experiment, 
the charged lepton emerging from the interaction vertex travels this distance within a time 
scale ts ~ 10~^ s, leading to an energy uncertainty above SE ~ h/ts ~ 10~^ eV. Taking as 
an example Sm^ ~ 10"^ eV^; U ^ Eg ^ 100 MeV, it follows that 6m^/U ~ lO^^^gy < 6E. 
Therefore, the detection or absorption of the charged lepton near the source, on a time scale 
much smaller than 1/A cannot resolve the energy difference between the mass eigenstates 



and the approximation (2.31) is justified 



Another approximation we can use in (2.26) is i7i ~ ^^2 — f2, for Q ^ A. Absorbing the 
phase factors e~*'^-'*s/2 j^to the definition of the states |z/j), the time evolved disentangled 
state is then approximately given by 
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(2.32) 
for ts < 27r/A. 

The state inside the brackets is identified as the usual quantum mechanical state that is 
time evolved from the "electron" neutrino state, which is prepared initially at ts- From 
this analysis, we see that there are two conditions required for the disentangled neu- 
trino state to be identified with the familiar quantum mechanical state: 5m?' /Vt -C 1 
and ts ^ tosc ~ l-n'^jbrn? . The former is always satisfied for neutrinos with 5vn? ~ 
10^'^ — 10~^eV^, fi > few MeV, while the latter is fulfilled if the charged lepton produced 
at the source is either observed or stopped at a beam-dump near the production region in 
long-baseline experiments. 

The latter condition implies that the neutrino state is disentangled before oscillations 
can occur. In a long-baseline experiment this is achieved if the charged lepton, which is 
entangled with the neutrino at the production vertex, is measured or stopped near the 
production region. Therefore, we conclude that an important condition for the familiar 
quantum mechanical description of oscillations to be reliable is that the quantum state must 
be disentangled on time scales much shorter than the oscillation time. 

It is clear from this discussion that the precise value of ts is irrelevant as long as ^5 ^ tosc- 
Experimentally, ts is of the order of the time of flight of the charged lepton in the production 
region at the source namely a few 10~* s. 

C. Transition amplitudes and event rates 



in, pit-ts) 



The number of charged lepton events with momentum Q produced at the source, at time 
ts is given by 



n,iQ,ts) = {^cits)\aliQ)a,m'^eits)) = {VeiQ,ts)\VciQ,ts))- 



(2.33) 



For ts ^ tosc = 27r/A and fi ^ A, using the approximations leading to (2.32), we obtain 
the charged lepton differential detection rate at the source 
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where at the last step we have replaced the diffraction function by the delta function. This 
can be justified as follows. For ts ~ 10~^ s, the width of this function (the resolution) 
in energy is ~ 10"'^ eV. Since the typical energy in a long-baseline experiment is > 40 — 
100 MeV, the error incurred in replacing the diffraction function by a delta function is smaller 
than one part in 10^^. 

We can also obtain the transition amplitude for the disentangled state to produce a final 
charged lepton and another W particle at the far detector at time to, where to — ts^L and 
L is the baseline. This is obtained from the time evolution of the disentangled state with 
the full Hamiltonian from the time of disentanglement ts until the time to, at which the 
detection measurement is carried out by overlapping the time evolved state with the final 
state. 

The transition amplitude is given by 



A, 



a->/3 



{W{kD),h{R 



D 



-iH(tD-ts) 



VeiQM) 



-,-iEDt 



" {W{kD),h{PD)\U{tD,ts)e'''''' \Ve{Q,ts)). 



(2.35) 



The disappearance and appearance amplitudes are then given by 
Ae^e = -/ n {2nf6{ks -Ps-ko- Pd) x 



cos^^ 



2n 



g "'l,P 2 



+ 



1,P 



sm 



sm 



{Es - n,^p) 



2n 



-in --^ 

_g "'2,P 2 



2,P 



{Es - n^^p)/2 
{Es 



sm 



{ED-n^p){tD-ts)/2 



sm 



"2,Py 2 



{Es - n,p)/2 



sm 



{Ed - f^i,p)/2 
{ED-n2p){tD-ts)/2 



{ED-n,p)/2 



(2.36) 



and 



A, 



e^fH 



-/ n {2iT)H{ks -Ps-ko- PD, 



sin 2^ 



X 



-in -^ 

g "'l,P 2 



2n 



1,-p 



sm 



-in --^ 

g "'2,P 2 



2Q 



2,P 



(^5-Mt 



{Es - l^i,p)/2 

[Es - ^2,p) ' 
{Es-n,p)/2 



sm 



sm 



{ED-n,p){tD-ts)/2 



sm 



{En - f^i,p)/2 

{ED-n2^p){tD-ts)/2 

{En - n,^p)/2 ^ 



(2.37) 
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with P = ks — ps, and 11 is given by 



16 V^ E^ E^ El El 

ks ko PS Vd 



(2.38) 



the labels D, S refer to (far) detector and source respectively. 



Implementing the same approximations leading to the factorized state (2.32), namely 



O 3> A and t^A ^ 1, we find the disappearance and appearance probabilities 

Ve^eito) = (^^y i2nf V 5iks - PS -kD~PD)27Tts6iEs- Up) 



cos^ e fUx, t, A) + sin^ 9 ft (x, t, A) + 2 cos' 6 sin' 6 cos(tA)/+(x, t, A)/_(a;, t, A) 



' e-^^iV' 



A'^D) 



7^0x2 



2n 



p 



{27if V 6{ks -ps-ko- Pd) 2vr ts 6{Es - Qj 



sin' 26 



/^(x,t,A) + /!(x,t,A)-2cos(tA)/+(x,t,A)/_(x,t,A)L 



(2.39) 



where t = to — ts and x = Eo — Vtp. Here, 



-^•S' ~ ^ks ~ PS ' ^D - i^t + -C/. 



kn+^PD^ 



and f± are given by Eq. 2.17 



(2.40) 



In the long time limit, using f±{x,t,A) — )• 2nt5{x ± A) their product is given by (2.19), 
and we find 



' e— s>el,'- 



D) 



9 



2n\2 



^^^^ (27rf V6{ks- PS -kD-PD)tsS{Es-Qp) 
cos^ et6{x + A) + sin^ et5{x-A) + 2 cos' 9 sin' 6 



Ve^,,{t 



D 



7'n\2 



2fi, 



(27r)^ V (5(A;5 - Ps - ko - Pd) ts SiEs - n 



sin(2tA) 1 
2A 2 

sin' 26 



6{x + A) + 6{x- A) 



t6{x + A)+tS{x-A)-2 



sin(2tA) 1 
2A 2 



P^ 4 
6{x + A) + (5(a: - A) 



(2.41) 



where P = k^ + pd. 

These transition probabilities feature the two time scales ts and t = to — ts associated 
with the measurements at the near and far detector. They also feature energy conserving 
delta functions associated with the different mass eigenstates. 

There is a further simplification when Q ^ A. In this regime, when the probabilities 



(2.41) are integrated over a smooth density of states, the delta functions corresponding to the 



mass eigenstates yield the density of states at values Er> = f2 =p A. In typical experiments, 
where ^2 ~ 100 MeV and 5m' ~ 10~^eV', the density of final states must vary dramatically 
near fl to resolve the small interval A, with A/Q < 10~^^. Therefore, understanding the 
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probabilities as being integrated with a smooth final density of states insensitive to the mass 
difference, we can approximate S{x±A) ~ S{x). In this case, we can approximate the above 
expressions by 



Ve^e{t 



D 



l^Il\^ 



2n, 



(27r)5 V 5{ks -ps~kn- Pd) ts 5{Es - np) 5{Ed - ^p) 



cos^ e + sin^ e 



5, ,,sin(2tA) 
+ 2cos2^sin2^ — ^— ^ 



(2.42) 



Ve^^^itD) 



7'n\2 



2fi, 



(27r)^ V 6{ks -ps-ko- Pd) ts 6{Es - Qp) 6{Ed - Qp, 



sin(2tA) 
2A 



sin^ 26 



(2.43) 



The product 6{Es — ^p) S{Ed — Qp) is an approximate energy conservation at both produc- 
tion and detection vertices, where we have neglected A, which is twice the energy difference 
between the mass eigenstates. 

Further insight can be gained by obtaining the phase space distribution of the number of 
charged leptons / = e, /i at the far detector 



(27r 



,3 dNr 
(Px d?PD 



ni{pD,tD) = {Ve{Q,tD,ts)\alipD)ai{pD)\Ve{Q,tD,ts)). (2.44) 



Here, 



\Ve{Q,tD,ts)) = e-'^°*^f/(ti5,ts)e'''°*^|Ve(Q,is)) 



(2.45) 



is the neutrino state disentangled at ts near the source and has been time-evolved until it is 
detected at the far detector at time tn- Not surprisingly, since we only keep terms that are 
up to order g"^, the time evolved state contains a single lepton Fock state. We find that 



{2ixf^^— = V,^,{to) ; {2nf ^ 



d^x d^p 



D 



d^x d^p 



' e-^fiV^ 



D 



e-^^J.y'^D), 



(2.46) 



with the probabilities Ve^eito) and Ve^^ito) are given by ( 2.42|2.43 ). 

Taking the time derivative with respect to t^, we obtain the differential charged lepton 
event rates at the far detector 



(2vr) 



dN. 



FD 



g^Ii^^ 



(27r)5 V 5{ks -ps-ko- Pd) ts 5{Es - np) 5{Ed - n, 



d^xdtd^PD V2i7p- 

cos^ e + sin^ e + 2 cos^ 9 sin^ 9 cos(2t A) 



(2.47) 



(27r^3 '^ 



d^xdtd^PD \2Vl 



g'liV 



(27r)S V 5iks -ps-ko- Pd) ts 5iEs - np) 5{Eu - n, 



1 - cos(2tA) 



sin^ 29 
2 

(2.48) 
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Remarkably, these rates can be simply written as 



(2vr) 



dN^ 



FD 



d?x dt d^pD 



(27r 



dN^ 



d'^x d^Ps 



a^P 



{t)dT 



up^-W Ipi 



(2.49) 



where we have used the expression (2.34) for the differential charged lepton event rate at 
the source and integrated in ts, 



kn PD t" 



^{ks -Ps-ko- Pd) S{Ed - fij 



(2.50) 



is the charged lepton production rate from the reaction i^^ — )■ Wlj3 for a flavor neutrino 
of energy ^2 and Pa^p{t) are the disappearance (a = /3) or appearance (a ^ (3) transition 
probabilities obtained from the usual quantum mechanical calculations of Rabi-oscillations, 



1 - sin2(20) sin^ 






P. 



e^ffj, 



sin^(2e) sin^ 



5rn^ 



(2.51) 



The remarkable aspect of the final result (2.50) is the factorization of the different pro- 



cesses contributing to the far detector event rate^ namely the number of events at the source 
multiplies the quantum mechanical transition probability which in turn multiplies the pro- 
duction rate at the vertex in the far detector. This factorization is a distinct consequence 
of the two time analysis, of the disentanglement of the neutrino near the production region 
along with the approximations invoked in the resolution of the energy conserving delta func- 
tions. We emphasize that the factorization in terms of the quantum mechanical transition 
probabilities Pa^^it) only applies to the detection rate defined by taking the time derivative. 
The total number of events per phase space volume also factorizes but not in terms of the 
quantum mechanical transition probabilities but in terms of their integral in time. 



D. Possible corrections for short baseline experiments 



This analysis also suggests that potentially important corrections may arise in short 
baseline experiments for much lower energy. In this case, it may occur that the time scale 
of disentanglement is of the same order as the oscillation scale Vt/5m? so that ts^rn"^/^ ~ 1 
and the sine functions cannot be factored out of the quantum state as in (2.32). In this case, 
the resulting disentangled state is 



\Ve{Q,t)) oc IsinO |z^2p) ^ 



AD. p{t-ts) 



/_ (x, ts, A) + cos I z/^ p) e 



-in p{t-ts) 



/+(x,ts,A) 



(2.52) 

which is not of the familiar quantum mechanical form of a coherent superposition of the 
mass eigenstates with the usual cos ^, sin ^ amplitudes but with the "diffraction functions" 
multiplying these factors. If the width of these functions is of order 211 /ts ^ 2A, interference 



term in the probability will feature the product (2.19) which would lead to an extra inter- 
ference terms of the form sin(t5A)/t5A multiplying the sin(2tA)/2A in (2.42|2.43). Such 
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factor yields an extra modulation with energy which may yield phenomenologically interest- 
ing modifications in the interpretation and analysis of data. Although the typical time scale 
ts for disentanglement near the source region is approximately the same for short and long 
baseline experiments, this discussion applies to the possibility in which ts is not too small 
compared to the oscillation time scale as could be the case in short-baseline experiments. For 
example at LSND the source region is a few meters across and the baseline is ~ 30 m, and 
at KARMEN where the baseline is ~ 17.5 m. The baseline and beam energy are designed so 
that at least an oscillation takes place along the baseline. However in these short-baseline 
experiments, the distance over which the neutrino state is disentangled from the charged 
lepton, a few meters, is not too much smaller than the baseline, consequently ts may not be 
much smaller than an oscillation time. 

The potential phenomenological importance of these effects merit their further study. 



E. Wave packet description 

As mentioned earlier in the paper, we restricted ourselves to plane-wave states, the main 
reason being to discuss the main concepts within the simplest setting. Although we postpone 
a thorough discussion of the more technical aspects including the wave packet description to 
further study, we can provide physical arguments that allow the extrapolation of the main 
results obtained above to the wave packet case. 

Consider that the initial particle, the charged leptons measured at the near and far detec- 
tors and the final particle are all described by wave packets with a macroscopic localization 
length scale a of the order of the typical scale of the detectors (a few meters) localized at the 
source and detectors respectively. These are sharply localized in momentum. The transition 



amplitudes between wave packets are now obtained from (2.35 2.37) by convolution with the 
wave functions of the initial state, the charged lepton measured near the source and those 
of the final state at the far detector. 



The disentangled neutrino state (2.26) is now described as a propagating wave packet 
with a typical spatial localization length also of 0{a), namely of macroscopic scale. These 
mildly localized wave packets lead to both momentum and energy uncertainty of the same 
order AE ~ Ap ~ h/a ~ h/ts ~ 10"''' eV. These uncertainties in energy and momentum are 
much larger than typical values for the neutrino energy differences 5m? /Vt ~ 10~"'^^eV. As 
this wave packet propagates, the mass eigenstates will slowly separate because they feature 
different group velocities, and coherence will be maintained as long as the separation between 
the wave packets is much smaller than their localization lengths (neglecting dispersion). For 
a baseline L and typical energies Vt ~ 100 MeV, the separation at the far detector is oc 
(A/fi) L ~ 10~^*^L. Therefore, coherence is always maintained for macroscopic localization 



lengths in terrestrial experiments. The detection amplitudes (2.36|2.37) and probabilities 



( 2.42|2.43 ) vanish until the time when the front of the packet overlaps with the far detector, 
namely for to — ts ^ L, with an uncertainty of 0{a), which for long-baseline experiments 
is always a <^ L. The total number of events at the far detector is found by integrating 



the rate (2.50) from t = L, which is when the front of the wave packet overlaps with the 
detector until t = L + a when the wave packet has completely passed through the detector. 
For L ^ a and neglecting the separation of wave packets the total number of events at the 
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far detector is obtained by replacing (2.49) by 



rfiVj^ _ dN^ 



d^x d^pn d^x d^p, [^"-^^^^^ ^] dV^^^wi, . (2.53) 

The delta functions will also be integrated with the wave packet profiles and volume factors 
are replaced by a^. Furthermore, at distances from the source ^ a, the disentangled neutrino 
wave packet propagates as a spherical wave. Therefore, for far-detectors in the "radiation 
zone," the probability will be suppressed by a geometric flux factor 1/L^ [2ll|25]. Obviously, 
there are uncertainties of 0{a) in this argument, but these are irrelevant for long-baselines 



L ^ a. The a multiplying the total number of events in (2.53) is expected: it is the "total 
interaction time" during which the wave packet interacts with the detector. Obviously, these 
factors modify the overall normalization. However, the distortion of the energy spectrum in 
disappearance and the appearance probabilities for long-baseline experiments are determined 
by Pa^j3{L, E) (for short baseline there are potential corrections as discussed above). While 
we deem these arguments describing the wave-packet scenario to be physically sound, we 
will provide a thorough technical discussion of this case based on the results obtained above 
in the plane wave approximation in future work. The analysis presented above for the 
dynamics of disentanglement and propagation provides the fundamental basis to include the 
wave packet description, which is obtained from simple convolution of these results with 
wave packet wave functions. 

III. CONCLUSIONS AND FURTHER DISCUSSION 

In appearance and disappearance experiments, neutrinos produced by a charged current 
vertex at the source are entangled with the charged lepton. We study the dynamics of this 
quantum state directly as a function of time. 

If the charged lepton (or daughter particle) produced at the source is not measured, 
tracing out this degree of freedom yields a density matrix for the neutrino. In the mass 
basis, the time evolution of the diagonal density matrix elements describe the production of 
mass eigenstates from the decay, whereas the off-diagonal density matrix elements exhibit 
the oscillations resulting from the interference of mass eigenstates and are a measure of 
coherence. We find that coherence remains up to a time of the order of the oscillation time 
scale during which diagonal and off diagonal matrix elements are of the same order. 

The measurement of the charged lepton near the source (or its stopping at a nearby 
"beam dump") disentangles the neutrino state, and it is the further time evolution of this 
disentangled state with the total Hamiltonian that leads to the production of charged lep- 
tons at the far detector. Thus, the process of production and detection in long-baseline 
experiments involves two different time scales: the measurement of the charged lepton near 
the source determines the first time scale at which the neutrino state is disentangled, while 
the measurement of the charged lepton at the far detector is the second time scale, which is 
much longer than the first one in long-baseline experiments. 

In this time dependent description, we establish that the usual quantum mechanical state 
emerges if the disentanglement of the charged lepton produced with the neutrino occurs on 
time scales much shorter than the oscillation scale ~ E/Sra^. Under these circumstances, 
the event rate at the far detector factorizes in terms of the usual quantum mechanical 
probability if the final density of states is insensitive to the difference in the energies of the 
mass eigenstates in the ultrarelativistic limit Sm'^/E. 

17 



Although in this article we focused on the study of the dynamics for the case of plane 
waves to exhibit the main results within the simplest setting, we provide physically moti- 
vated arguments to extrapolate the results to the case of wave packets, analyzing the effect 
of localization of the initial and final state wave functions at the near and far detector on 
the detection process. Combining the results obtained for plane waves with a wave packet 
description yields the total number of events detected in terms of the usual transition prob- 
ability as a function of baseline. A deeper and more detailed analysis of wave packets and 
localization aspects will be provided in forthcoming work. 

The analysis presented here also suggests that there could arise potentially interesting cor- 
rections in the case of short baseline experiments such as MiniBoone and/or LSND with base- 
lines of ~ 500, ~ 30 meters and typical energies ~ IGeV; ~ 30MeV, respectively, wherein 
the disentanglement of the neutrino could occur on time scales of the same order as the 
oscillation scale. In this case, the disentangled state will differ from the usual quantum me- 
chanical coherent superposition and this difference introduces extra modulation with energy 
and could be important in the final detection rate. Such possibility could be relevant in the 
interpretation and analysis of data, in particular masses and mixing angles and certainly 
merits further detailed study. 
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